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Abstract 

Colloidal systems (and protein solutions) are often characterized by 
attractive interactions whose range is much smaller than the particle size. 

When this is the case and the interaction is spherical, systems obey a gen- 
eralized law of correspondent states (GLCS), first proposed by Noro and 
Frenkel [ J.Chem.Phys. 113, 2941 (2000) ]. The thermodynamic prop- 
erties become insensitive to the details of the potential, depending only 
on the value of the second virial coefficient B2 and the density p. The 
GLCS does not gcncrically hold for the case of non-spherical potentials. In 
this Letter wc suggest that when particles interact via short-ranged small- 
angular amplitude patchy interactions (so that the condition of only one 
bond per patch is fulfilled) it is still possible to generalize the GLCS close 
to the liquid-gas critical point. 
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In colloidal and in protein systems, the interaction potential is often short- 
ranged, i.e. small as compared to the particle size. The small-range interaction 
brings some peculiar behaviour both on the system thermodynamics and on the 
dynamics. For example, the gas-liquid phase separation becomes metastable 
with respect to crystallization [Tj (and hence a proper equilibrium liquid phase 
is missing). In addition, for very short ranges, the kinetic arrest (glass) line 
becomes reentrant and two different glass phases appear |21 E] . The presence of 
short-range attractions is also invoked to explain the so-called "crystallization 
slot" in the phase diagram of globular proteins |4]. 

A common property unifies all spherically symmetric short ranged attractive 
potentials, independently from their actual shape. Indeed, Noro and Frenkel [5] 
showed that the thermodynamical properties of systems interacting with short- 
range attractive potentials are all equivalent if scaled by the proper variables, 
i.e. they obey a generalized law of correspondent stated (GLCS). They showed 
that the virial coefficient can be used as scaling variable for the strength of the 
interaction. Recently, the GLCS has been shown to arise from the fact that, 
due to the short-range of the interaction, each interacting pair of particles (a 
bond) contributes independently and equally to the partition function [S] . 

The above considerations are valid for centrosymmetric potentials and can 
not be straightforwardly extended to non-spherical cases [71 [S], i.e. when interac- 
tions are patchy and strongly directional. In the case of molecular fluids, patchy 
interactions are relevant in network- forming systems, like silica [SJ [TD] and wa- 
ter [m [T2] and in all associating fluids [13 E] in which the hydrogen bond plays 
an important role. For colloidal systems, this interest is justified by recent ad- 
vances in the synthesis and characterization of patchy particles [HI \W\ [18] 
or functionalized colloidal particles [TH]. For what concern proteins, it is weU 
established that the interactions are intrinsically directional [201 [HI [22 and that 
their description in term of isotropic potentials represents only a first-order ap- 
proximation. 

Recent progress in understanding the thermodynamic of patchy colloidal parti- 
cles has been based on application of the thermodynamical perturbation theory 
developed by Wertheim 23J . The theory, which does not account for the geom- 
etry of the patches, assumes that each patch acts as an independent interacting 
unit. In addition, the theory neglects the possibility of close-loops of bonds. 
Despite these approximations, several predictions of the theory have been nu- 
merically confirmed [T31 IMl [21] . It has been shown that (for the case of patchy 
interactions) the number of independent interacting patches (the valence) is the 
key ingredient in controlling the phase diagram of the system. On lowering the 
valence, the liquid-gas critical point shifts to smaller and smaller densities, so 
that liquid states of vanishing density (empty liquids) become accessible [24] . 
The Wertheim expression for the bonding free-energy is a function only of the 
bond probability and of the valence. In this respect, the theory suggests that 
systems with the same valence should behave similarly, if the bond probability 
is the same. Since the bond probability is related to the chemical bonding con- 
stant [13j |25| , and since for large attraction strengths the chemical constant is 
proportional to the second virial coefficient [20], the Wertheim theory suggests 
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that universality based on the virial scaling can be recovered also in the case of 
patchy interactions when valence is preserved. 

To address the issue of a possible generalization of the Noro-Frenkel scaling 
we have extended the investigation concerning the location of the critical point 
for the patchy model introduced by Kern and Frenkel [7j, for several values of 
the attraction range and of the width of the patchy attractive region. The choice 
of this potential is motivated by the fact that angular and radial properties can 
be independently modified. More precisely, the two body potential is defined 
as: 

^.(r,,) = u'^{r,,)f{{n,,}) (1) 

where u"'^ (vij ) is an isotropic square well term of depth mq and attractive range 
(7 + A and f{{^ij}) is a function that depends on the orientation of the two 
interacting particles {^ij}- The diameter of the particles and the depth of the 
square well has been chosen as units of length and energy respectively, i.e. a — 1 
and uq = 1. Each particle is characterized by M identical patches. A patch 
a is defined as the intersection of the surface of the sphere with a cone with 
half-opening angle that has the vertex in the center of the particle and it has 
the axis directed toward the direction u„. The angular function /({ilij}) is 
defined as 
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where Tij is the direction of the vector that joins the centers of the two inter- 
acting particles and a {(3) some patch belonging to the particle i (j) . In practice 
two particles interact attractively if, when they are within the attractive dis- 
tance CT -I- A, two patches are properly facing each other. When this is the case, 
the two particles are considered bonded. Decreasing A reduces the range of the 
attraction whereas reducing 9 diminishes the angular size of the patches. In 
the limit A the model goes toward the patchy Baxter limit. In the limit 
cos 9 1, the patch goes to the point limit. 

In the present work we focus on M=3, M=4 and M=5 patches, located on the 
surface of the particle as shown in the cartoon of Fig. [T] Differently from previ- 
ous studies [7j, we consider values of A and 9 such that, due to steric reasons, 
each patch is involved simultaneously in only one pair interaction, i.e. those val- 
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ues fuUfiUing the condition sin6' > 2{a'+A)^ ■ Under this single-bond per patch 
condition, the number of patches coincides also with the maximum number of 
possible bond per particle. 

The Kern-Frenkel potential posses an analytical expression for the second 
virial coefficient: 

^ = l-X^((l + A)3-l)(eV^-l) (3) 
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where x = A/ ^~^°^(^) jg the percentage of surface covered by the attractive 
patches and the temperature is measured in reduced units, i.e. ks = 1. Here 
is the hard-sphere component of the virial coefficient. To calculate the 
location of the gas-liquid critical point we perform grand canonical Monte Carlo 
(GCMC) simulations [26], complemented with histogram reweighting techniques 
to match the distribution of the order parameter p—se with the known functional 
dependence expected at the Ising universality class critical point J^7\. Here e 
is the potential energy density, p the number density and s is the mixing field 
parameter. We did not performed a finite size study, since we are only interested 
in the trends with the range A and the angular size of the patches 9. We have 
studied systems of size L = 6 for M = 4 and 5 and L = 7 for M = 3, where 
L is the side length of the cubic simulation box. For each studied M — using 
the methods described in |25] — we calculated the critical temperature Tc and 
density pc for values of A between 0.119 and 0.01 (at fixed cos 61 = 0.92) and for 
cos9 between 0.895 and 0.99 (at fixed A = 0.119). The results are summarized 
in Table [H 

We start by analyzing the results as a function of A, at fixed cos6' = 0.92. 
Fig. |2] shows that Tc decreases with A while pc increases. The A dependence of 
Pc can be conveniently described by Pc(A) = Pc{0) /(l-t-0.5A)^, a functional form 
which suggests that pc would be constant if measured using as unit of length the 
average distance between two bonded particles (1 -I- 0.5A). The resulting pc{0) 
extrapolated value provides an estimate of the corresponding patchy Baxter 
model Pc- Fig. [2]-(a) shows also that the Tc dependences are apparently well 
described by iso-i32 lines. Each M is characterized by a different B2 value, 
enforcing the existence of a GLCS for each valence. As a confirmation, we 
evaluate the values of 82/62^ at the critical point for the patchy particle model 
studied in Ref. [24]. The resulting values for M — 3,4 and 5 are respectively 
Ba/Sf ^ = -28.12, -4.95 and -2.78, very similar to the values reported in the 
present study. We also note that the B2 differences between different M values 
are much larger than the variation with A at constant M. Hence at a zero-th 
order approximation, when the repulsive part of the potential is complemented 
by localized patchy interactions, B2 can be considered as a scaling variable of a 
GLCS in the single-bond per patch condition. 

A closer look to the actual B2 values (see Table [T|) shows that a small trend 
in the B2 values is present, which is hidden in the logarithmic transformation 
relating B2 to T (see Eq.|3|. Still, the B2 differences between different M values 
are much larger than the variation with A at constant M. This suggests that, 
at a zero-th order approximation, B2 can be considered as a scaling variable 
of a GLCS in the single-bond per patch condition. Thus, provided that the 
geometry of the patches is such that their number coincides with the maximum 
valence, B2 carries the information on the valence of the patchy interaction 
potential. These results suggest that, statistically, configurations with the same 
Boltzmann weight are generated under an isotropic scaling (to change the inter- 
particle distances preserving the same bonding pattern) and a simultaneous 
change of both A and T such that B2 remains constant. It is interesting to 
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observe that the GLCS for the isotropic square well is fulfilled for values of the 
range smaller then 0.05 HH], whereas for the three models discussed here 
GLCS appears to hold for longer ranges. 

Fig. [3] shows Tc and Pc, this time at A = 0.119, as a function of cos9. 
For values of cos 9 larger than the one reported in the figure, crystallization 
is observed within the simulation time, preventing the possibility of evaluating 
the critical parameters. This effect suggests that the liquid-gas critical point is 
metastable, as in the case of spherical short-range potentials. The observation 
of crystallization inform us already that on reducing 9, some bonding patterns 
acquires a larger statistical weight. Indeed, differently from the previous case, 
small deviations from a GLCS are observed in T (see lines Fig. [3]-(a)) and 
not only in B!^- We attribute these deviation to the fact that by changing 
the angular part of the potential, the statistical relevance of specific bonding 
patterns varies. In other words, it is not possible to vary the orientation of 
the particles to preserve the bonding on changing 9, i.e. it is not possible 
to perform the operation equivalent to rescaling the distances to preserve the 
bonding pattern on changing A. We expect that the breaking of the scaling 
will be enhanced at state points far from the critical region where an extensive 
bonding pattern is present, i.e. low T and large p. Concerning the 9 dependence 
of Pc, we note that it increases with cos 6', for the case M = 5 and M = 4 while 
it weakly decreases for M — 3. We have no clear arguments for interpreting the 
Pc{d) trends, except for the fact that the observed dependence (more significant 
for M = 4 and M = 5), suggests a non-trivial coupling between the angular 
correlation induced by bonding and the density. For these two M values, small 
bonding angles appear to require smaller average distances between next-nearest 
neighbor particles implying larger densities. 

The argument behind the quasi-validity of the GLCS for each M class are 
based on the fact that the relevant role is played by the bonding pattern [5], 
which is supposed to be statistically identical for all members of the class along 
corresponding states. Hence, the number of bonds at the critical point should 
be similar. To double check this statement we also report in Table [T] the bond 
probability pg, defined as the potential energy at the critical point normalized by 
the energy of the fully bonded system. Such quantity is indeed constant for each 
M value, in agreement (and strongly supporting) the possibility of defining a 
different GLCS for each valence class. We also note that, within each M class, 
the variation of are smaller than the one of B2 , suggesting that the bond 
probability may result in a better scaling variable than B2. A small trend in 

is only observed in the 9 dependence. We also note that this observation is 
in agreement with the Wertheim theory [SS] and with the identification of the 
bond free energy as the appropriate scaling variable for the spherical case [6]. 

In summary, we have provided evidence that different short-ranged non- 
spherical potentials, but with the same number of single-bond patches, essen- 
tially obey a GLCS. The condition of a single-bond per patch requires that both 
the attraction range and the angular size of the patches are small. Breakdowns 
of the GLCS can be expected for potentials which differ in their angular part, 
especially for very small angular sizes (cos 9 — *■ 1), since — in conditions of 



5 



extensive bonding — the statistical weight of closed loops of bonds becomes 
significantly affected by the angular patch size. 
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M 


A 


cos 9 


Tc 


Pc 


Pi 




3 


0.050 


0.92 


0.1076 


0.151 


0.727 


-23.67 


3 


0.070 


0.92 


0.1113 


0.150 


0.732 


-24.92 


3 


0.119 


0.92 


0.1180 


0.134 


0.721 


-26.71 


3 


0.119 


0.895 


0.1252 


0.136 


0.737 


-28.29 


3 


0.119 


0.92 


0.1180 


0.134 


0.721 


-26.71 


3 


0.119 


0.95 


0.1069 


0.131 


0.726 


-25.13 


3 


0.119 


0.97 


0.0968 


0.129 


0.726 


-23.92 


4 


0.010 


0.92 


0.1140 


0.326 


0.656 


-4.02 


4 


0.050 


0.92 


0.1392 


0.306 


0.652 


-4.32 


4 


0.070 


0.92 


0.1468 


0.306 


0.641 


-4.24 


4 


0.090 


0.92 


0.1511 


0.289 


0.642 


-4.65 


4 


0.119 


0.92 


0.1573 


0.276 


0.644 


-4.92 


4 


0.119 


0.895 


0.1682 


0.267 


0.635 


-5.74 


4 


0.119 


0.92 


0.1573 


0.276 


0.644 


-4.92 


4 


0.119 


0.94 


0.1484 


0.307 


0.667 


-3.85 


5 


0.010 


0.92 


0.1259 


0.393 


0.567 


-2.41 


5 


0.050 


0.92 


0.1570 


0.387 


0.588 


-2.67 


5 


0.070 


0.92 


0.1653 


0.374 


0.585 


-2.81 


5 


0.090 


0.92 


0.1723 


0.362 


0.581 


-2.90 


5 


0.119 


0.92 


0.1808 


0.348 


0.577 


-3.03 


5 


0.119 


0.895 


0.1959 


0.333 


0.570 


-3.53 


5 


0.119 


0.92 


0.1808 


0.348 


0.577 


-3.03 


5 


0.119 


0.95 


0.1601 


0.372 


0.589 


-2.22 



Table 1: Critical point properties for all studic^d mod(Js, labeled by the values 
M, A and cos 9. The different columns indicate the critical temperature Tc, 
critical density pc , the bond probability and the reduced value of the second 
virial coefficient at the critical point. The estimated errors for each 

of these quantities are: ±0.0005(rc); ±0.007(pc); ±0.005(pg). The error in B2 
arises from the error in Tc and differs for each point due to the non-linear relation 
between B2 and T. The field-mixing parameter s is always smaller than 0.08. 
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Figure 1: Pictorial representation of the three different patch geometries con- 
sidered in this work. The soUd angle of each patch is 27r(l — cos^), where 6 is 
the cone semi-angle. 
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Figure 2: Critical temperature Tc (a) and critical density Pc (b) as a function 
of the range A, at fixed patch angular size cos9 = 0.92. Lines in (a) correspond 
to constant values of B2: specifically B2/B^^ = -3.03 for M = 5, B2/B^^ = 
-4.92 for M = 4 and Sa/Sf = -26.71 for M = 3. Lines in (b) are pc{A) = 
Pc(0)/(1 + O.SA)-^, where Pc{0) (shaded point) is a fitting parameter. The fit 
values are pdO) = 0.41 for M = 5, pc(0) = 0.33 for M = 4 and pc{0) = 0.16 for 
M = 3. Open symbols are simulation results. 
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Figure 3: Critic;al temperature Tc (a) and critical density pc (b) as a function of 
the angular patch size cos 6*, at fixed range A = 0.119. Lines in (a) correspond 
to constant values of B2: specifically -Bz/Sf = -3.53 for M = 5, B2/B^^ = 
-5.74 for M = 4 and B^/B^^ = -28.29 for M = 3. Open symbols are 
simulation results. 
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